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Kinetics of collision processes with linear mixing rules are investigated analytically. The velocity 
distribution becomes self-similar in the long time limit and the similarity functions have algebraic or 
stretched exponential tails. The characteristic exponents are roots of transcendental equations and 
vary continuously with the mixing parameters. In the presence of conservation laws, the velocity 
distributions become universal. 
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Collision processes underlie fundamental phenomena 
such as heat transport in gases Q and mixing in fluid 
flows In closed systems, conservation of mass and mo- 
mentum implies Maxwellian velocity statistics How- 
ever, there are physical systems such as granular media 
Q and atomic collisions 5] where energy or even momen- 
tum 1^ are not conserved. The reason may be that only 
a subset of the system is considered, that not all degrees 
of freedom are measured (2D imaging of 3D systems), 
or that acoustic or other excitations are ignored. Non- 
Maxwellian velocity statistics are found in granular gases 
, colloids , turbulent flows |^ , and laser cooling . 

Random collision processes are a framework for study- 
ing the role of conservation laws, and demonstrate how 
anomalous velocit y st atistics emerge when conservation 
laws are relaxed [H El El El. Motivated by this, 
we consider binary collision processes with arbitrary lin- 
ear collision rules. While in the long-time limit velocity 
distributions are generically self-similar, there is a wide 
spectrum of possible behaviors. The velocity distribu- 
tions are characterized by algebraic or stretched exponen- 
tial tails and the corresponding exponents depend sensi- 
tively on the collision parameters. Interestingly, when 
there is energy or momentum conservation, the behavior 
is universal. 

Let us consider the most general linear collision law: 
when a particle of velocity ui collides with a particle of 
velocity U2, its post-collision velocity vi is given by 



vi = pui + qu2, 



(1) 



with p and q the mixing parameters. Special cases include 
elastic collisions (p = 0, q = 1), inelastic collisions {p + 
q = 1), the granules model {p + q < 1) |6|, the Kac Model 
(p^ 4- q^ = 1) 15], inelastic Lorenz gas {p = 0, q < I) 
|l^. and addition {p ^ q = 1) [l7l |. 

Further, we consider perfectly random dynamics: two 
randomly chosen particles collide according to The 
normalized velocity distribution P{v, t) obeys 



d_ 

dt 



P{v,t) = 



duidu2P(ui, t) P{U2, t) 



(2) 



X [5{v - pui - qu2) ~ 5{v - ui)] 



This Boltzmann equation is termed the Maxwell model in 
kinetic theory (the constant collision rate is set to unity) 
E3- The quadratic integrand in Eq. Q reflects the bi- 
nary nature of the collision process and the gain term 
reflects the collision rule ^ . The number density is con- 
served by Eq. / dvP{v,t) = 1. 

The convolution structure of the evolution equation 
suggests the Fourier transform F{k,t) = J dve^'^'" P{v,t). 
This quantity obeys the nonlinear and nonlocal equation 



d_ 

dt 



F{k, t) + F{k, t) = F{pk, t)F{qk, t). 



(3) 



This closed equation is amenable to analytical treat- 
ment. Moments of the velocity distribution, Mn{t) 



J dvv"P{v, t), obey a closed hierarchy of equations 



dt""- 



n-l 

E 

m— 1 



(4) 



with the shorthand notation A„ = 1 — p" — (7" . These 
equations are solved recursively with Mo{t) = 1. 

We are interested in the long time limit and we seek 
similarity solutions of the form 



P{v,t) 



*$(we° 



as t 



00 . 



(5) 



Equivalently, the Fourier transform has the similarity 
form F{k,t) /(fce^"*). This function satisfies 



azr{z) + f{z)^f{pz)f{qz). 



(6) 



The similarity function may include both a regular and 
a singular component f{z) — fi-cg{z) + /sing(-z) with 

/rog(z) = J2n -^^Z"- Normalization sets /q = 1. 
The leading small-z behavior of the singular component, 
fsing{z) ~ z'^, reflects an algebraic tail of the velocity 
distribution 



<i>(u)) 



w 



(7) 



as w 



-> 00. Substituting the leading singular behavior 
1 ^ z" into the governing equation JHJl yields a 
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FIG. 1: The phase diagram for type-I Scahng. Shown are 
equi-!^ contours as a function of the mixing parameters. 

relation between the scaling parameter a and the mixing 
parameters p and q 

a = v-^K. (8) 

There are two types of scaling solutions, depending on the 
average initial velocity Afi(O): type-I scaling (Afi(O) 7^ 0) 
and type-II scaling (Mi(0) — 0). 

Type-I Scaling. The first moment varies exponentially 
with time, Mi — e"^^* , according to the rate equation 
+ XiMi = (its initial value can be set to unity). 
The small wave-number behavior of the regular com- 
ponent of the Fourier transform is therefore F{k,t) ^ 
1 -|- ike~^^^. When ly > 1, this component dominates 
over the singular component. Therefore, f{z) ^ 1 + iz 
and the scaling parameter is a = Ai. When < 1, the 
singular component dominates over the regular one, so 
the parameter a is not obvious. There is a spectrum 
of possible a's depending on ly according to the "disper- 
sion" curve ||SJ). This curve has a maximum a = -^X^ a,t 
v given by v-^Xi, = X^. We argue that this maximum is 
actually realized by the dynamics. Intuitively, the selec- 
tion of the extremum point maximizes the typical wave 
number oc e"'. The scaling parameter is therefore 

/ p^lni + fl'^lni i/<l; 
y \ — p — q i^>l. 

The exponent v characterizing the algebraic tails is the 
root of the transcendental equation (jHJ with a given by 
@ . Explicitly, the equations are p" In ^ + q'^ In = 1 
and 1 — p'^ — q'^ — v{l — p — q) ior V < 1 and v > 1, 
respectively (Fig. 1). 

Algebraic tails exist as long as the exponent v is fi- 
nite. The exponent diverges, — > 00, in the limiting 
cases p-\- q = 1, p = 1, and q = 1, defining the triangular 



FIG. 2: The phase diagram in case II. 

region S (Fig. 1). In this domain, the singular compo- 
nent disappears and the scaling function f{z) is analytic. 
Moreover, the velocity distribution ^{w) has sharp tails 
and all of its moments are finite. Outside the region S, 
the exponent u is always finite. It varies continuously 
as a function of the mixing parameters and it vanishes, 
v —^ 0, when p or q 0. 

The curve pln^-|-gln| = I, marking the case v = I, 
separates two kinds of behavior. When 1/ > 1 the first 
moment characterizes the velocity distribution. In the 
complementary case, the typical velocity docs not follow 
from the (integer) moment behavior. This dichotomy is 
reminiscent of similarity solutions of the first and second 
kind p5j . Interestingly, extremum selection determines 
the typical velocity and thereby the velocity distribu- 
tion when < I. Extremum selection similarly governs 
the speed and the shape of nonlinear waves in reaction- 
diffusion problems [21|. Indeed, in terms of the variable 
Inv, the similarity solution Q is nothing but a travelling 
wave P{lnv,t) <&(lnti -I- at). 

Type-II scaling. When the average initial velocity van- 
ishes, the initial variance can be set to unity. From the 
moment equations the variance varies exponentially 
with time, M2{t) — e^'*'^*. The small- wave number be- 
havior F{k,t) = 1 — ifc^e"^^* dominates when v > 2 
and consequently, the second moment characterizes the 
velocity distribution, a — ^A2. Otherwise, the singular 
component governs the behavior as above and the scaling 
parameter is 

r n"^ In i + fl'^ In i v < 2: , , 

\\{\-p^-q^) v>2. ^ ' 

The exponent v characterizing the algebraic 
tails is the root of the transcendental equa- 
tion In + In -17 = 1 for v < 2 and 
1 -p"^ -g"^ = 1(1 -p2 _q2) for ,y > 2 (Fig. 2). It 
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FIG. 3: Self-similarity in type-I scaling. Shown is $(it)) versus 
w for the case p = q — A {a — .2). In the inset, the tail is 
compared with the theoretical prediction i/ = 4.88636. 
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FIG. 4: Self-similarity in type-II scaling. Shown is ver- 
sus w for the case p = q — .6 {a = 0.14). The inset compares 
the tail with the theoretical prediction u — 6.66937. 



diverges in the vicinity of the curves p = 1, q = 1, 



and p 



1. Inside this region S, the similarity 



solution f{z) is regular and the velocity distribution 
^{w) has sharp tails. Generally, similarity solutions are 
symmetric, as f{z) = f{—z). 

The diverging exponent v indicates that the large ve- 
locity tail is a stretched exponential, rather than alge- 
braic, in the region S. Indeed, the Fourier transform 
/(z) exp(— z^) for large z is compatible with the gov- 
erning equation ij^l when 



= 0. 



(11) 



This in turn suggests a stretched exponential behavior 



~ exp (- 



-w ' 



(12) 



with 7 = — ^ for w — > oo. As was the case for the 
exponent j/, /i is a root of a transcendental equation and 
consequently, the exponent 7 > 1 varies continuously 
with the mixing parameters. Its minimal value 7 = 1 is 
attained along the region boundaries p — 1 and q — 1. 
The tail is Gaussian, 7 = 2, along the type-II boundary 
A2 = and the exponent diverges, 7 ^ 00, along the 
type-I boundary Ai = 0. 

Self-similarity holds whether the typical velocity 
shrinks or grows with time, i.e., regardless whether a 
is positive or negative. For type-I scaling the velocities 
shrink (grow) with time when Ai > (Ai < 0) and simi- 
larly for type-II scaling. 

Asymptotically, sufficiently small moments of the ve- 
locity distribution are governed by the typical velocity. 
Otherwise, the moment behavior follows from the hierar- 
chy of evolution equations Q 



exp(— noft), n < v; 
exp(— A„t), n> v. 



(13) 



Sufficiently small moments of the velocity distribution 
exhibit ordinary scaling behavior while sufficiently large 
moments exhibit multiscaling asymptotic behavior 23|. 
By multiscaling, we refer to moment ratios such as 
Mn/Mj'^^ that diverge asymptotically. 

Numerical simulations confirm the theoretical findings. 
In the simulations, randomly chosen pairs of particles 
undergo the collision process The number of parti- 
cles was 10^ and the velocity distributions were obtained 
from an average over 10 independent realizations. Sim- 
ulation results corresponding to flat initial distributions 
with support in [0 : 1] and [—1 : 1] are shown in Fig. 3 
and Fig. 4, respectively. 

There are a number of special cases worth highlighting. 

1. Addition (p = q — 1): This integrable case nicely 
demonstrates how similarity solutions emerge. The 
transformation G{k,t) = 1/F{k,t) reduces the (local) 
Ricatti Eq. @ into -§^G{k,t) + G{k,t) = 1. The Fourier 

transform reads F{k,t) = {1 + [i^T^fc) - l]e*}"^ with 
Fo{k) = F{k,t = 0). Indeed, the small wave number 
behavior of the initial distribution Fo{k) dictates the 
asymptotic behavior. When Mi(0) = 1, type-I scaling 
occurs, with the similarity solution f{z) = [l — iz]~^ and 
<^{w) = e^"" for u; > 0. When Mi(0) = 0, type-II scahng 
occurs, with the similarity solution /(z) = [1 + ^z^]~^ 

and ^{w) — exp(— V2|w|). One can verify that these 
similarity solutions satisfy ^ with a — I and 1/2 for 
type-I and type-II, respectively. 

2. Kac Model (p^ + q^ = 1 ): For type-II scaling, energy is 
conserved since A2 = 0. In this case, the velocity distri- 
bution approaches a steady state, a = 0. The equation 
(O has the solution f{z) = exp(— z^/2) and the velocity 
distribution is Maxwellian <I>(u)) — (27r)^-'^/^ exp(— w^/2). 
For type-I scaling, energy is not conserved and the veloc- 
ity distribution is no longer universal. However, it still 
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exhibits a Gaussian tail. 

3. Inelastic Maxwell Model (p + q = I): For inelastic 
collisions, the total momentum is conserved, Ai = 0. Us- 
ing the Galilean transformation v v — Afi(O), the ini- 
tial momentum can be set to zero and so the behavior 
is always of type-II. The exponent = 3 is the root 
of the equation = ■|A2. In this particular case, an 
explicit solution can be found f{z) = (1 -I- |z|)e~l^l or 
^{w) = ^{l + 'w^)~^ Interestingly, when there is a 
conservation law (either momentum or energy), the sim- 
ilarity solution is independent of the mixing parameters. 
4- Lorentz gas (p = 0): This case corresponds to inelastic 
collisions with massive scatterers. The evolution equa- 
tion is linear, ^P{v,t) +P{v,t) = q'^ P{vq~^ ,t). It is 
useful to consider the stochastic process the velocity un- 
dergoes V vq vq^ —>•••. The number of collisions 
is distributed according to a Poisson distribution with 
mean equal to time t. Therefore, the velocity distribu- 
tion is [23 

''(M)--f:S;^^'o(;^). (14) 

In other words, the variable Inw is Poisson distributed 
with mean equal to t In q, so a finite number of stan- 
dard deviations away from the mean \nv is Gaus- 
sian distributed. Thus, the tail of the distribution is 
log-normal, P{v) ~ exp[— (lnw)^/(2i Ing)]. In the limit 
1/ "> 0, no similarity solutions emerge and all mo- 
ments of the velocity distribution exhibit multiscaling, 
Mn{t) = M„(0) exp(— A„t). We conclude that the Hnear 
collision process is an effective mixing mechanism. 
No matter how small either of the mixing parameter is, 
eventually, the binary collision process alters the nature 
of the velocity distribution. In other words, nonlinearity 
provides the mechanism for the self-similar behavior. 



The usual physical particle collision region is as- 
sociated with 0<p-fg<l(fi| and restitution 
coefficient < q — p < 1. Other applications may 
involve values of p and q outside of this region. We 
implicitly assumed that the mixing parameters are 
positive (j), 9 > 0) but the behavior easily extends to 
the other three quadrants in the p — q plane. Consider 
the first two moments Mi(i) = Mi(0) exp(— Ait) and 
M2{t) = [M2(0) - cMi2(0)]e-^2* + cMf{0)e-^^^' with 
c = {2pq)/ {X2 — 2X1). The first moment governs the 
second moment (M2 ~ M^) when A2 > 2Ai, i.e., in the 
circular domain {p — 1)^ + {q — 1)^ < 1, that is entirely 
contained in the first quadrant. Thus, type-I scaling 
occurs only in the first quadrant. Type-II scaling occurs 
in the other three quadrants, regardless of Mi(0). 

We tacitly assumed that all moments are finite ini- 
tially. Consider initial distributions with a leading small- 
k behavior of the type 1 — FQ{k) ~ k'^", competing with 
/sing(z) ~ z'^ . The initial conditions govern the asymp- 
totic behavior when i^q < ly and thus, a — ly^^X^f^. This 
generalization of the previous results applies for both 
type-I scaling (i/q — 1) and type-II scaling {i^o — 2). 

In closing, random and linear mixing results in self- 
similar velocity distributions. Nonlinearity is responsible 
for this scaling and extremum selection may govern the 
behavior. The velocity distributions have either algebraic 
or exponential tails, with nontrivial characteristic expo- 
nents. Every possible algebraic tail and every faster than 
exponential decay constitutes the spectrum of behaviors. 
Conservation laws play a crucial role, as the velocity dis- 
tribution becomes universal when physical quantities (ei- 
ther energy or momentum) are conserved. 
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